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Abstract - This paper addresses the detection of low velocity 
target from a shipborne HFSWR (High Frequency Surface 
Waves Radar). In a shipborne configuration, as a 
consequence of the ship motion first order of Bragg lines is 
spread in the useful Doppler interval where detections of 
ship target are expected to appear. Space time techniques 
like STAP, usually studied for airborne radar purposes are 
candidate to reduce the effect of ship motion. Influence of 
pitch and roll components is also discussed. Simulations of 
STAP techniques are given, based on a conventional 
architecture of processing. Estimation of the clutter 
covariance matrix is studied with respect of the 
characteristics of the radar waveform. Interleaving between 
the classical HFSWR waveform and a wide band learning 
waveform is proposed to enable the covariance estimation. 


I. INTRODUCTION 

Surface waves propagation is well known for its 
benefit to radar detection of target beyond the horizon. 
There are a lot of references and published papers on 
applications of HF surface waves radar (HFSWR) to 
oceanography (sea current measurement) , target 
detection and tracking. Most of the HFSWR systems are 
onshore radar and provide interesting detection 
performances against low altitude targets like missile, 
stealth aircraft or vessels. However, when the radar has 
to be mounted on a shipborne platform, additional 
problems have to be solved, such as HF antennas 
integration, E.M. compatibility and motion 
compensation. In a shipborne configuration the first 
order Bragg lines is spread, making the detection of ships 
more difficult than in a onshore configuration, because 
the Doppler frequency of these targets are expected to 
appear in the spreading clutter domain. This effect is as 
more problematic as velocities of the ship (hosting the 
radar) and the ship target present the same order of 
magnitude. 

Shipborne HFSWR configuration is similar to 
airborne early warning configuration where low altitude 
target echoes are embedded in the ground spread clutter. 
To counter this effect, signal processing techniques like 
DPCA (display phase center antenna) and - more 
recently - STAP (Space Time Adaptive Processing) can 
be used to restore the performances. During the last ten 
years, the abundance of studies and works to optimise 


STAP architecture was coming from airborne radar 
domain, and purpose of this paper is study whether these 
algorithms can be adapted and finally implemented on a 
shipborne HFSWR to help the detection of slow targets. 

Purpose of the first part of this paper is to remind 
the basement of physics of Bragg scattering. In a second 
part, simulations of Bragg clutter in a motion configuration 
are presented. The third part describes the possible 
architecture of a STAP processor with respect to the 
HFSWR specificity. Simulations are presented to assess 
the feasibility of STAP and the problem of clutter 
covariance matrix is discussed. 


II. SEA CLUTTER IN HFSWR 

Sea clutter consists of 1 st order and second order of Bragg 
lines, within an interval of a few Hz. The 1 st order, as 
received by an onshore radar, is distributed on 2 resonant 
frequencies : 

F B ragg+/- = -2 UO.VA +/- 0.1 (F HH z) 1/2 

with : 

Uo : radar steering vector , 
v : current velocity, 

F MHz ’ c/X =w/2n (X = wave length) 

The first order reflectivity is equal to : 
a a \w) = a{S(co+ 2 k 0 .v - 1 /gj2^f).r(6>) + S(co+ 2 k 0 .v + ^g|2£ 0 | ).r(6> +7i)} 

where the magnitude coefficient a depends on the sea state 
and the frequency. For a given carrier frequency (for 
example 15 MHz), this coefficient reaches an asymptotic 
level of - 23dBm 2 /m 2 when the wind velocity exceeds a 
certain threshold (for example 3 m/s @15 MHz). 

Calculation of second order of Bragg lines is based on a 
dual series of wave numbers [k,kl] presenting a resonant 
combination according to the radar wavelength. A 
simplified computation technique [1] derived from Barrick 
theory can be used : 


0-7803-7871-7/03/$17.00 © 2003 IEEE 


504 


Radar 2003 



Report Documentation Page 

Form Approved 

OMB No. 0704-0188 

Public reporting burden for the collection of information is estimated to average 1 hour per response, including the time for reviewing instructions, searching existing data sources, gathering and 
maintaining the data needed, and completing and reviewing the collection of information. Send comments regarding this burden estimate or any other aspect of this collection of information, 
including suggestions for reducing this burden, to Washington Headquarters Services, Directorate for Information Operations and Reports, 1215 Jefferson Davis Highway, Suite 1204, Arlington 

VA 22202-4302. Respondents should be aware that notwithstanding any other provision of law, no person shall be subject to a penalty for failing to comply with a collection of information if it 
does not display a currently valid OMB control number. 

1. REPORT DATE 

14 APR 2005 

2. REPORT TYPE 

N/A 

3. DATES COVERED 

4. TITLE AND SUBTITLE 




5a. CONTRACT NUMBER 

Use of STAP techniques to enhance the detection of slow targets in 

5b. GRANT NUMBER 

Miipuui lie nr o vv xv 



5c. PROGRAM ELEMENT NUMBER 

6. AUTHOR(S) 




5d. PROJECT NUMBER 





5e. TASK NUMBER 





5f. WORK UNIT NUMBER 

7. PERFORMING ORGANIZATION NAME(S) AND ADDRESS(ES) 

Office National dEtudes et Recherches Aerospatiales (ONERA) Chemin 
de la Huniere et des Joncherettes, F-91761 Palaiseau 

8. PERFORMING ORGANIZATION 

REPORT NUMBER 

9. SPONSORING/MONITORING AGENCY NAME(S) AND ADDRESS(ES) 


10. SPONSOR/MONITOR’S ACRONYM(S) 





11. SPONSOR/MONITOR’S REPORT 
NUMBER(S) 

12. DISTRIBUTION/AVAILABILITY STATEMENT 

Approved for public release, distribution unlimited 

13. SUPPLEMENTARY NOTES 

See also ADM001798, Proceedings of the International Conference on Radar (RADAR 2003) Held in 
Adelaide, Australia on 3-5 September 2003., The original document contains color images. 

14. ABSTRACT 

15. SUBJECT TERMS 

16. SECURITY CLASSIFICATION OF: 

17. LIMITATION OF 
ABSTRACT 

uu 

18. NUMBER 
OF PAGES 

6 

19a. NAME OF 
RESPONSIBLE PERSON 

a. REPORT 

unclassified 

b. ABSTRACT 

unclassified 

c. THIS PAGE 

unclassified 


Standard Form 298 (Rev. 8-98) 

Prescribed by ANSI Std Z39-18 






ship 


(j (2 \co) = \C\ 2 4{2K) 5 k*\d 2 k X T.S(mk)S(pkl) {Tik,^)] 2 

m=± 1 p =±1 

. S(Aa>+2k 0 V-mjg\k\ -p^fg\k^\) 

with 


k).{k Q .k x )- 

where S(.) represents the sea wave gravity spectrum. 
Second order of Bragg line also depends on frequency 
and sea state. In the upper HF band (15-30 MHz) the 
level of the second order is much more sensitive to sea 
state than in lower band (below 10 MHz), where sea 
surface can be considered as "smooth according to the 
wavelength". Level of second order can be from -60 to - 
34 dBm 2 /m 2 (for a coherent integration time of 100 s) 
depending on the wavelength and the seastate [2]. 

Then, in a motionless configuration the Bragg 
spectrum presents the shape shown on figure 1. Spectrum 
has been computed for a coherent integration time of 90 
seconds, at 15 MHz, for a wind of 8m/s. The Pierson 
Moscovitz gravity spectrum was considered. Targets of 
interest (ships) are endo-clutter with a Doppler between - 
1Hz and 1Hz (radial velocity lower than lOm/s @ 15 
MHz). 


2 (k 1,)} 


and 


|C| 2 =l/(2;r) 4 


T(k,k\) = 


{(*o 



Fig. 1. : Bragg spectrum, 15 MHz, wind 8 m/s; 
no current; wind vector parallel to radar axis 



- 70dB -30 dB OdB 


Fig. 2 . Example of HFSWR detection at 16 MHz , coastal radar 


III. CLUTTER IN SHIPBORNE HFSWR 

Motion consists of translation, pitch and roll, as shown on 
figure 3. We consider the following set of parameters (for 
a sea state 5) relative to a large ship (frigate class) 
equipped with a dynamic stabilisation system (transfer of 
ballast): 

Roll : magnitude +/- 6.5°, period 12 seconds 
pitch : magnitude +/-10 0 , period 6 seconds 
velocity (V) : lOm/s 

Due to linear motion, sea echoes in each azimuth 0 suffer 
from an additional Doppler shift equal to: 

Af=2V/A, .sin(6) 

Thus, clutter is distributed along an oblique ridge in the 
plane [Af, sin(0)]. 

Effect of pitch and roll has to be computed. A dual rotation 
matrix is considered and the additional phase distortion 
calculated between each antenna and any current position 
M (figure 3). 


In a motionless configuration (for example coastal 
radar), ship echoes are detected through Doppler 
processing and two cases are possible: 

If the target Doppler shift coincides with one of the 
Bragg lines, carrier frequency has to be changed to 
provide another chance to detect the target. Else, 
waiting for Doppler changing (i.e. several minutes 
for ships) allows to detect the target. 

If the target detection is disturbed by the second 
order of Bragg clutter, detection enhancement can 
be obtained from improving the range resolution, as 
far as a wider frequency bandwidth can be available 
for radar detection purpose. 

Figure 2 shows an experimental result obtained from a 
coastal HFSWR system at 16 MHz. Target is a 
commercial ship (oil tanker class). 


M (scatter) 



The map of figure 4 gives the Doppler response of any 
point M (of polar coordinates [0 ,r]) in a wide angular 
sector from -90° to 90° with reference to the normal axis 
of the ship. Only phase distortion has been considered. 
Pitch and roll create additional parallel clutter ridge of 
lower magnitude The configuration is equivalent to have a 
translation motion with a modulated velocity. Effect is 
more significant for high value of squint angle 0. 
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Doppler - Hz 

Fig. 4. Azimuth - Doppler response 

From the Bragg spectrum and the Doppler response 
studied previously it is now possible to compute the 
effect of motion on a global landscape of sea surface 
scatters. Assuming no correlation between the statistics 
of signals coming from different directions we can write: 

S(f) = \D(6)Rdop(ej)SBra g8 (f).dd 


where D( 0) is the directional pattern, S bragg (f) is the Bragg 
spectrum(l st and 2 nd orders) and R dop (0,f) the Doppler 
response computed above. 

Figure 5 gives the result at 15 MHz and 25 MHz, for a 
sea state equal to 5. Several curves are superimposed in 
order to understand the contribution of each component of 
the motion to the global distortion. Thus, with reference to 
the onshore configuration, translation creates a spread of 
the 1 st order Bragg line. Pitch and roll without translation 
create spurious at a level below the 2 nd order. When the 
ship is moving, pitch & roll effect does not modify 
significantly the shape of the clutter curves (see the black 
and red dotted lines). 

All these curves of course take into account a 
weighting window (Hanning) to reduces the Doppler 
processing secondary lobes and avoid a wrong 
interpretation of clutter spurious. 

This first analysis shows that in our case (small 
pitch & roll deviations) the main effect is due to 
translation. Then, priority will be given to algorithms such 
as STAP, matched to linear motion and able to restore the 
detection of slow targets originally embedded in the spread 
sea clutter. 




Fig. 5 . Simulation of motion effect on Bragg spectrum, including pitch and roll 
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IV. 


USE OF STAP TECHNIQUES 
A. Introduction 

STAP consists of adaptive techniques which allow 
to reduce the influence of clutter spread spectrum. 
These techniques are used in airborne domain to 
enhance the detection of small target , also to 
achieve the detection of ground target (GMTI). 
Compared to airborne domain, HFSWR present 
several differences : 

frequency band and waveform (PRF, duty 
cycle, unambiguous in range and velocity) 
clutter phenomenology (Bragg scattering) 
velocity of the platform, 

In a conventional STAP approach, a space time 
weighting vector W is searched as a solution of a 
minimisation problem : 

L(W,A,) = W h RW+ X ,(T h W- 1) , where W,X are 
unknown and R is the space time covariance matrix 
The well known solution presents the following 
form : 

W = R-'t(t h R-'T f 1 y = - R 5 

opt \ / r-pH r>l'T"' 


where X is the data vector, T is steering vector of 
the target (for a linear array, the phase is equal to 
27i// 1 [i a sin(6 target )a + i r V target ], where a is the 
array step and Tr the repetition period, [i a ,i r ] being 
the index of antennas and recurrences) 

We can notice that R^X operates a filtering of 
input signals with respect of the clutter covariance 
matrix. Digital beamforming and Doppler 
processing are derived from multiplication by J H . 
In fact, this last operation can also be done more 
easily through a 2D-FFT and the normalisation 
term TR d T can be omitted in a simplified 
approach, if a CFAR detection is applied after 
STAP filtering. 

Figure 6 shows the conventional architecture of a 
STAP processor. The receiving antennas are 
gathered in subarrays (in our case for example, 16 
antennas in 4 subarrays of 4 elements). 


The sub array architecture is required to reduce the 
complexity and achieve the accurate estimation of the 
covariance matrix R. Feasibility of this estimation will be 
discussed in a further section. Data from each sub array are 
delivered at different times, following uniform taps x, 2 t, 
..Mx. The delay x represents a multiple of the recurrence 
period, for example 500 ms in the case of HFSWR. This 
delay is corresponding to a first coherent Doppler 
integration of 500 pulses, assuming a classical PRF of 
lKHz (figure 7). 


Sub array 1 



x = 500 ms, M=9, N= 4 or 8 


Fig. 6 . STAP conventional architecture 


Received signals on one channel 



-0HH3—-Q-D-0- 


Post-processing (20 snapshot = 90s) 


Fig. 7 . Time sequences for STAP processing 
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B. STAP Simulation 


Simulation of STAP has been done in the following 
conditions : 

carrier frequency : 15 MHz 

radar signals are considered in a Doppler 

interval of -1 Hz a +1 Hz following a first 

coherent integration of 500 ms duration(figure 

7). 

Bragg lines are considered with the following 
arbitrary magnitude : 65% for negative line, 
25% for the positive line, and 10% for the 
centre line (zero Doppler); for each line, a 
random gaussian process has been generated to 
represent the sea scatters statistic. Not any 
correlation is assumed to exist between the 
negative and positive line. These two 
assumptions (gaussian and no correlation) 

should be verified through experimental 

investigations. 

4 or 8 antennas channels have been considered 
from a receive passive antenna installed on the 
broadside of the ship. The antenna is 120 
meters length . 

9 taps (of 500ms)are considered. 

size of the covariance matrix is [72 x 72] or [36 

x 36] depending on the number of sub arrays. 

STAP filter is applied to signals along a duration of 
9 x 500 ms = 4.5 sec. 

A post processing of 20 STAP results is made, so 
that the total processing time is equal to 90 seconds. 
The ship velocity is equal to 10 m/s. 

Simulation results are shown on figure 8. 
Target presents a Doppler shift of 0.2 Hz and an 
azimuth angle of -17°. Target amplitude has been 
adjusted to present - in a exo-clutter assumption - a 
signal to thermal noise ratio of 13 dB. 

As the filter is not normalised by l^R' 1 ! , 
we can see in the region of clutter the presence of 
residual echoes (false alarms). However the side 
lobes of the main clutter region has been reduced 
making possible the detection of the target. 

Post-processing has been considered in 
this simulation as non coherent. Coherent 
processing is possible but requires to compensate 
the pitch and roll from one snapshot of to the 
following (every 4.5s). This compensation can be 
done using the pitch & roll measurements provided 
by a position and orientation system. Commercial 
on the shelves equipment provide angular accuracy 
better than 1° and can be used successfully. A 
coherent integration (compared to a non coherent 
one) allows to improve the signal to noise ratio of 
3dB and to get a better Doppler discrimination. 



dop Hz 



dop Hz 

Fig. 8. Simulation of STAP technique 


C. Estimation of Covariance matrix 

In the simulation shown above, covariance matrix has been 
estimated around the target cell, without any constraint on 
availability of range samples. Estimation of covariance 
matrix is subjected to gather enough samples in range to 
make this estimation accurate. Usually, estimation of a 
covariance matrix of size [NM,NM], 2NM samples are 
needed to approach the optimal performances (-3dB loss), 
where N is the number of taps and M the number of space 
channels. Assuming 9 taps (delays) and 4 space channels, 
72 range cells are needed. Considering a range resolution 
of a few kilometres to ten kilometres (a common value for 
HF radar), the range interval required for estimation should 
be very large, up to 200 or 300 km. In this interval, 
variation of sea current and sea state can make risky and 
inaccurate the covariance estimation. Then, the solution 
consists of using an wider band (thinner range resolution) 
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waveform dedicated to covariance estimation and a 
conventional waveform for radar detection. These 
two waveforms are interleaved as shown on figure 
9. Repetition period is doubled. The wide band 
waveform is a sparse pulse waveform where the 
pulse duration is the same as for the original 
waveform; frequency is changing from pulse to 
pulse according to a set of frequencies previously 
selected with respect of occupied channels. 
Covariance matrix can be estimated in frequency 
domain more easily than in time domain. Assuming 
a range interval of 15 km for covariance estimation, 
pulse duration is equal to lOOqs. For a PRF of 
500Hz, up to 250 frequencies can be selected. For a 
covariance matrix size of [36,36], 72 channels are 
needed for the estimation. Assuming that 50% of 
the channels are not clear, the total required 
bandwidth (B cov , for covariance estimation) is equal 
to 720 kHz around the centre radar frequency. For a 
radar frequency of 15 MHz, estimation has to be 
performed between 14.64 MHz and 15.36 MHz. 
Within this bandwidth, dispersion due to variation 
of Bragg lines (from 0.38262 Hz to 0.39192 Hz ~ 9 
mHz) is negligible for a STAP applied on a 
duration of 4.5 seconds. Antenna calibration needs 
also to be considered in this wide bandwidth; as for 
Bragg lines, strong dispersion of antenna response 
is theoretically unexpected if the length of the array 
(-100 m) is smaller than cl B cov , and small 
variations can be compensated through calibration. 


Radar 

Detection 

fo — 


100 JLXS 


Coy. Radar Coy. 

Estimation Detection Estimation 


F] 

// 

fo 


Fk 


// 




2 ms 

- ->j 




-~-m --> 

500 ms 


radar 



<-> 

B cov =Few 100 kHz 


Fig. 9. HFSWR waveform matched to STAP 


Due to the high reflectivity of clutter 
(compared to level of desired targets), the average 


power per frequency channel, for the 'learning' waveform 
can be very low, thus compensating the inconvenient of 
such a wide band transmission in the HF band. 

Other solutions, to decrease the complexity of 
STAP implementation are on study. In particular, 
emergence of sub-optimal techniques presents a great 
interest [3]. 


CONCLUSION 

This paper addressed the effect of motion on HFSWR and 
proposed the use of STAP techniques to improve the 
detection of slow targets. 

Effect of linear motion, which produces a spread 
of Doppler echoes, can be compensated by use of STAP 
techniques in the same way as in airborne radar domain for 
GMTI purpose. STAP is also robust to small magnitude of 
pitch and roll, as considered in this paper. However, the 
case of higher angular deviation could be problematic in so 
far as linear techniques (like covariance filtering) cannot 
compensate rotation effects. 

STAP has been discussed and illustrated from 
simulation results, assuming a gaussian clutter distribution 
and no correlation between the Bragg lines. Next steps of 
work must address experimentation with collection of real 
data in order to refine the properties of the sea clutter (in 
particular to address the correlation aspect of Bragg lines) 
and assess the performances of STAP algorithms. 

Feasibility of STAP implementation has been 
discussed, taking into account the specificity of HFSWR 
waveform. To counter the rarefaction of range cells (due to 
the narrow band available on this kind of radar) we 
proposed to estimate the covariance matrix in frequency 
domain, from a set of sparse frequency available in the 
neighbourhood of the radar carrier frequency. This 
technique might be also verified and validated from real 
data in a next experimental step. 
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